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ABSTRACT 

The giant planet atmospheres exhibit alternating prograde (eastward) and retrograde (westward) jets of different speeds and widths, 
with an equatorial jet that is prograde on Jupiter and Saturn and retrograde on Uranus and Neptune. The jets are variously thought to 
be driven by differential radiative heating of the upper atmosphere or by intrinsic heat fluxes emanating from the deep interior. But 
existing models cannot account for the different flow configurations on the giant planets in an energetically consistent manner. Here a 
three-dimensional general circulation model is used to show that the different flow configurations can be reproduced by mechanisms 
universal across the giant planets if differences in their radiative heating and intrinsic heat fluxes are taken into account. Whether the 
equatorial jet is prograde or retrograde depends on whether the deep intrinsic heat fluxes are strong enough that convection penetrates 
into the upper troposphere and generates strong equatorial Rossby waves there. Prograde equatorial jets result if convective Rossby 
wave generation is strong and low-latitude angular momentum flux divergence owing to baroclinic eddies generated off the equator is 
sufficiently weak (Jupiter and Saturn). Retrograde equatorial jets result if either convective Rossby wave generation is weak or absent 
(Uranus) or low-latitude angular momentum flux divergence owing to baroclinic eddies is sufficiently strong (Neptune). The different 
speeds and widths of the off-equatorial jets depend, among other factors, on the differential radiative heating of the atmosphere and 
the altitude of the jets, which are vertically sheared. The simulations have closed energy and angular momentum balances that are 
consistent with observations of the giant planets. They exhibit temperature structures closely resembling those observed, and make 
predictions about as-yet unobserved aspects of flow and temperature structures. 



1. Introduction 

Among the most striking features of the giant planets are 
the alternating zonal jets. As shown in Fig. 1, Jupiter 
and Saturn have prograde equatorial jets (superrotation) 
that peak at m and --200-400 m s~\ depend- 
ing on the vertical level considered. Uranus and Neptune 
have retrograde equatorial jets (subrotation) that peak 
at ~100ms~^ and ^150-400ms"^. Jupiter and Sat- 
urn have multiple off-equatorial jets in each hemisphere; 
Uranus and Neptune have only a single off-equatorial jet 
in each hemisphere. Despite decades of study with a va- 
riety of flow models, it has remained obscure how these 
different flow configurations come about (Vasavada and 
Showman 2005). 

Existing models posit as the driver of the flow ei- 
ther the differential radiative heating of the upper at- 
mosphere (e.g., Williams 1979, 2003b) or the intrinsic 
heat fluxes emanating from the deep interior (e.g., Bus se 
1976; Heimpel et al. 2005; Aurnou et al. 2007; Chan and 
Mayr 2008; Kaspi et al. 2009). However, none of these 
models can account for the existence of equatorial super- 
rotation on Jupiter and Saturn and equatorial subrotation 
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on Uranus and Neptune with radiative heating, intrin- 
sic heat fluxes, and other physical parameters consistent 
with observations. 

For example, deep-flow models that posit intrinsic 
heat fluxes as the sole driver of the flow can generate 
equatorial superrotation, but they use heat fluxes more 
than 10^ times larger than those observed (e.g., Heimpel 
and Aurnou 2007). They generate equatorial subrotation 
only with intrinsic heat fluxes even stronger than those 
for which they generate superrotation (Aurnou et al. 
2007), although the intrinsic heat fluxes on the subrotat- 
ing planets (Uranus and Neptune) are weaker than those 
on the superrotating planets (Jupiter and Saturn). The 
relevance of such deep-flow models is further called into 
question by the eddy angular momentum fluxes they im- 
ply. Their meridional eddy fluxes of angular momentum 
per unit volume (taking density variations into account) 
have a barotropic structure: they extend roughly along 
cylinders concentric with the planet's spin axis over the 
entire depth of the fluid, typically to pressures of or- 
der 10^ bar (e.g., Kaspi et al. 2009, their Fig. 10). But 
the eddy angular momentum fluxes inferred from track- 
ing cloud features in Jupiter's and Saturn's upper tropo- 
spheres indicate that the mean conversion rate from eddy 
to mean-flow kinetic energy is of order 10 ~^ W 
(IngersoU et al. 1981; Salyk et al. 2006; Del Genio et al. 
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Figure 1 Mean zonal velocities in the upper atmosphere of the giant planets from observations and simulations. Jupiter: ob- 
servations from the Cassini spacecraft (Porco et al. 2003) (orange line), and in simulation at 0.75 bar (dark blue line). Saturn: 
observations from the Voyager spacecraft (orange line), from the Hubble Space Telescope (HST) (Sanchez-Lavega et al. 2003) 
(green crosses), from the Cassini spacecraft at ^0.06 bar (magenta circles) and at ^0.7 bar (light blue squares) (Sanchez-Lavega 
et al. 2007), and in simulation at 0.1 bar (dark blue line). Uranus: observations from the Voyager spacecraft (orange circles), HST 
(orange crosses) (Hammel et al. 2001), the Keck telescope (orange squares) (Hammel et al. 2005), and in simulation at 25.0 mbar 
(dark blue line). Neptune: observations from the Voyager spacecraft (orange circles) and from HST (orange crosses) (Sromovsky 
et al. 2001), and in simulation at 25.0 mbar (dark blue line). Differences between the statistically identical northern and southern 
hemispheres in the simulations are indicative of the sampling variability of the averages. 



2007). If the observed upper-tropospheric eddy fluxes 
of angular momentum per unit volume extended unabat- 
edly over a layer of 50 km thickness (e.g., from about 
0.3 to 2.5 bar pressure on Jupiter, or from about 0.3 to 
0.9 bar pressure on Saturn), and if vertical zonal flow 
variations over this layer are not dramatic, the total en- 
ergy conversion rate would be about 0.5 Wm~^. This 
is already ~4% of the total energy uptake of the atmo- 
sphere from intrinsic heat fluxes and absorption of solar 
radiation for Jupiter, or '^11% for Saturn. But the limited 
thermodynamic efficiency of atmospheres allows only a 
fraction of the total atmospheric energy uptake to be used 
to generate eddy kinetic energy (Lorenz 1955; Peixoto 



and Oort 1992). The observations of Jupiter and Saturn 
therefore imply that eddy angular momentum fluxes can- 
not extend unabatedly over great depths and must have 
a baroclinic structure. Barotropic eddy angular momen- 
tum fluxes that extend to depths of order 10^ bar, with 
upper-atmospheric fluxes of similar scale and magnitude 
as those observed, are only possible in deep-flow models 
if the driving heat fluxes are several orders of magnitude 
greater than observed. 

Similarly, shallow-flow models that posit differential 
radiative heating as the driver of the flow can generate 
equatorial superrotation, but they require artifices such 
as additional equatorial heat or wave sources that have 
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no clear physical interpretation (e.g., Williams 2003b, a; 
Yamazaki et al. 2005; Lian and Showman 2008). It is 
unclear in those and other shallow-flow models (e.g., 
Scott and Polvani 2008) what physical characteristics 
distinguish the superrotating planets from the subrotat- 
ing planets.^ 

In Schneider and Liu (2009) (SL09 hereafter), we pos- 
tulated that prograde equatorial jets on the giant planets 
occur when intrinsic heat fluxes are strong enough that 
Rossby waves generated convectively in the equatorial 
region transport angular momentum toward the equator. 
Multiple off-equatorial jets, by contrast, form as a result 
of baroclinic instability owing to the differential radia- 
tive heating of the upper atmosphere. We introduced a 
general circulation model (GCM) and demonstrated with 
it that the postulated mechanisms can account qualita- 
tively for large-scale flow structures observed on Jupiter. 
Here we use simulations with essentially the same GCM, 
with closed energy and angular momentum balances that 
are consistent with observations, to demonstrate univer- 
sal formation mechanisms of jets on all the giant plan- 
ets. We show that the different flow configurations on 
the giant planets can be explained through consideration 
of the different roles played by intrinsic heat fluxes and 
solar radiation in generating atmospheric waves and in- 
stabilities. 

Section 2 briefly describes the GCM. Section 3 shows 
simulation results for Jupiter, Saturn, Uranus, and Nep- 
tune. Section 4 discusses the formation mechanisms of 
the jets in the simulations and confirms the postulated 
mechanisms through control simulations. Section 5 dis- 
cusses what the upper-atmospheric fluid dynamics, on 
which we focus, imply about flows at greater depth on 
the giant planets. Section 6 summarizes the conclusions 
and their relevance for available and possible future ob- 
servations. 



^Lian and Showman (2010) claim that different rates of latent heat 
release in phase changes of water may be responsible for superrota- 
tion on Jupiter and Saturn and subrotation on Uranus and Neptune. 
However, they impose latent heat fluxes at the lower boundary of their 
model that are not consistent with the observed energetics of the plan- 
ets. Similar to the simulations of Aurnou et al. (2007), they require 
stronger energy (latent heat) fluxes to generate subrotation than to gen- 
erate superrotation. For example, the latent heat fluxes are of order 
10-20 Wm~^ in their Jupiter and Saturn simulations and of order 
1500 W in their Uranus/Neptune simulation (Y. Lian, pers. com- 
munication, 2010). The latter are several orders of magnitude larger 
than the observed intrinsic heat fluxes or absorbed radiative fluxes (Ta- 
ble 1), which would have to drive any latent heat fluxes (energy would 
be required to evaporate the condensate that falls from the upper atmo- 
sphere into deeper layers). 



2. General circulation model 

With current computational resources, it is not feasible to 
simulate flows deep in giant planet atmospheres, where 
radiative relaxation times are measured in centuries and 
millennia, while at the same time resolving the energy- 
containing eddies in the upper atmospheres. Therefore, 
we focus on flows in the upper atmospheres, using a 
GCM that solves the hydrostatic primitive equations for 
a dry ideal-gas atmosphere in a thin spherical shell. The 
model is essentially that introduced for Jupiter in SL09, 
but here we use it to also simulate Saturn, Uranus, and 
Neptune.^ Parameters such as the planetary rotation rate, 
gravitational acceleration, and material properties of the 
atmosphere in each simulation are those of the planet be- 
ing simulated. The resolution in each simulation (T85 to 
T213 spectral resolution in the horizontal and 30 or 40 
levels in the vertical) is sufficient to resolve baroclinic 
instability and the energy-containing eddies in the upper 
atmosphere. The GCM and the simulations are described 
in detail in the appendix, and Table 1 lists the parameters; 
here we only give a brief overview. 

The GCM domain is a thin but three-dimensional 
spherical shell, which extends from the top of the atmo- 
sphere to an artificial lower boundary. The mean pres- 
sure at the lower boundary is 3 bar in all our simulations, 
to minimize differences in arbitrary parameters among 
them. Insolation is imposed as perpetual equinox with 
no diurnal cycle at the top of the atmosphere. Absorp- 
tion and scattering of solar radiation and absorption and 
emission of thermal radiation are represented in an ideal- 
ized way that is consistent with observations where they 
are available (primarily for Jupiter). Where radiative and 
other parameters are not well constrained by observa- 
tions or by knowledge of physical properties of the plan- 
ets, we set them to be equal to the parameters for Jupiter, 
again to minimize differences in unconstrained param- 
eters among the simulations. A dry convection scheme 
relaxes temperature profiles in statically unstable layers 
toward a convective profile with dry-adiabatic lapse rate, 
without transporting momentum in the vertical (see the 
appendix for details and for a discussion of this idealiza- 
tion). 

At the lower boundary of the GCM, a temporally 
constant and spatially uniform intrinsic heat flux is im- 
posed, with magnitude equal to the observed intrinsic 
heat fluxes (5.70, 2.01, and 0.04, and 0.43 Wm-^ for 
Jupiter, Saturn, Uranus, and Neptune). (The heat flux in 
the Uranus simulation corresponds to an observational 
upper bound.) Linear (Rayleigh) drag retards the flow 

^The Jupiter simulation here differs slightly from that in SL09 in 
that poorly constrained drag parameters in it are chosen to be the same 
as in the simulations of the other giant planets presented here. 
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Table 1 Simulation Parameters 



Parameter, symbol 


Jupiter 


Saturn 


Uranus 


Neptune 


Planetary radius, a (10^ m) 


69. 86 


c ^ ooa 

57.32'^ 


25.27 


o /I cr cr b 

24.55 


Planetary angular velocity, Q (10~^ s) 


1.7587 


1.5388 


1.U124 


1 rkoo /I b 

l.Uoo4 


Gravitational acceleration, g (ms~^) 


ZO.U 


1 n rb 




1 1 ob 
ii.Z 


Specific gas constant, R (Jkg~^ K~^) 


ooUo.oo 


/I ni /lb 


Q 1 /I n ob 


Q1 OT Tb 

oiy / . / 


/vuidDdLic exponent, rc 


9 /7 


9 /7 


9 /7 
Z/ ( 


9 /7 
Z/ / 


Specific heat capacity, Cp = R/k, (10^Jkg~^K~^) 


i.zo 


1 /1 1 
i.4i 


I.IU 


1 1 o 

i.iz 


Solar constant Fq (W m~^) 


oU. / 


1 /I oc 


Q 71 C 
O. / 1 


i.OZ 


Intrinsic heat flux (W m~^) 


0.7 


Z.Ul 


U.U42 


n /I Qoe 


JjUIIU. dlUCUU, / QQ 






u.ou 


n 9Qt> 
u.zy 


Single-scattering albedo, Co 


0.8 


0.8 


0.8 


0.8 


Solar optical depth at 3 bar, Tso 


3.0 


3.0 


3.0 


3.0 


Thermal optical depth at 3 bar, r/o 


80.0 


120.0 


60.0 


40.0 


Drag coefficient, ko (days"^ = (86400 s-^)) 


1/100 


1/100 


1/100 


1/100 


No-drag latitude, 00 


33° 


33° 


33° 


33° 


Horizontal spectral resolution 


T213 


T213 


T85 


T85 


Vertical levels 


30 


30 


40 


40 


Cut-off wavenumber for subgrid-scale dissipation 


100 


100 


40 


40 



"Guillot (1999); ^Lodders and B. Fegley (1998); ^Levine et al. (1977); ^Gierasch et al. (2000); 
^Guillot (2005); -^Hanel et al. (1981); ^Hanel et al. (1983) 



away from but not near the equator — a thin- shell repre- 
sentation of a magnetohydrodynamic (MHD) drag that 
acts at great depth (at pressures > 10^ bar), where the 
atmosphere becomes electrically conducting (Liu et al. 
2008). Drag in a deep atmosphere affects the angular 
momentum balance averaged over cylinders concentric 
with the planet's spin axis, so there is no effective drag 
on the flow in the upper atmosphere near the equator, 
in the region in which the cylinders do not intersect the 
layer of MHD drag at depth (SL09). Absent detailed 
knowledge of where and how the MHD drag acts and 
to rule out that differences among the simulations are 
caused by differences in the drag formulation, we chose 
the equatorial no-drag region to extend to = 33° lati- 
tude and the drag coefficient outside this region to be the 
same in all simulations. Section 4d discusses the effect 
of this drag formulation on our simulation results, and 
section 5 and the appendix provide further justification 
for it. 



We show simulation results from statistically steady 
states, which were reached after long spin-up periods; 
see the appendix for details. The northern and south- 
ern hemispheres in the simulations are statistically iden- 
tical, so differences between the hemispheres in figures 
showing long-term averages are indicative of the sam- 
pling variability of the averages. 



3. Simulation results 

a. Upper-atmospheric zonal flow 

Figure 1 shows the simulated mean zonal velocities near 
the levels at which cloud features from which the ob- 
served flows are inferred are suspected to occur: in the 
Jupiter simulation at 0.75 bar, corresponding to the layer 
of ammonia ice clouds on the actual planet (Atreya et al. 
1999); in the Saturn simulation at 0.1 bar, in a layer 
of tropospheric (e.g., ammonia) hazes (Sanchez-Lavega 
et al. 2007); in the Uranus and Neptune simulations at 
25 mbar, near the top of the stratospheric layers in which 
hydrocarbons would condense and form hazes (Gibbard 
et al. 2003). 

The simulations reproduce large-scale features of the 
observed flows in the upper atmosphere (Fig. 1). The 
Jupiter and Saturn simulations exhibit equatorial super- 
rotation, the Uranus and Neptune simulations equatorial 
subrotation. The equatorial jet in the Jupiter simulation 
has similar strength (^150 ms~^) and width as the ob- 
served jet. The equatorial jet in the Saturn simulation 
is stronger (^230 ms~^) and slightly wider than that in 
the Jupiter simulation, but it is weaker and narrower than 
the observed jet at a corresponding level on Saturn. The 
Jupiter and Saturn simulations exhibit alternating off- 
equatorial jets; they are broader than the observed jets 
but of similar strength. Especially in the Saturn simula- 
tion, the retrograde jets (except for the first retrograde jet 
off the equator) are broad and weak with speeds less than 
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10 m . They are more manifest as local minima of the 
zonal velocity than as actual retrograde jets. The Uranus 
and Neptune simulations exhibit a single off-equatorial 
jet in each hemisphere. The overall structure of the jets 
in the Uranus and Neptune simulations is roughly consis- 
tent with observations, but the equatorial jet in the Nep- 
tune simulation at the level shown ('^— 40ms~^)is con- 
siderably weaker than that observed. (However, the jet 
is stronger at higher levels in the simulation; see Fig. 5 
below.) 

In general, the prograde jets (or zonal velocity max- 
ima) are sharper than the retrograde jets (or zonal ve- 
locity minima), consistent with the zonal velocity max- 
ima being barotropically more stable (Rhines 1994). In- 
deed, meridional gradients of both absolute vorticity and 
quasigeostrophic potential vorticity at the levels at which 
the zonal velocities are shown in Fig. 1 are small near 
zonal velocity minima and are reversed near some of 
them (Fig. 2). (The quasigeostrophic potential vortic- 
ity is not shown for the Jupiter simulation because the 
stratification at the corresponding level is nearly stati- 
cally neutral, so that potential vorticity is not well de- 
fined; see Fig. 5 below.) The changing magnitude of 
the vorticity gradients between zonal velocity maxima 
and minima gives rise to a staircase pattern of absolute 
vorticity and potential vorticity as a function of latitude 
(Mclntyre 1982; Dritschel and Mclntyre 2008). Ab- 
solute vorticity gradients can reach about —2/3, partic- 
ularly in higher latitudes in the flanks of the minima; 
quasigeostrophic potential vorticity gradients are also re- 
versed near some of the zonal velocity minima, particu- 
larly in the Uranus and Neptune simulations, but they do 
not reach as strongly negative values as the absolute vor- 
ticity gradients. These features are roughly consistent 
with observations of Jupiter and Saturn (Ingersoll et al. 
1981; Read et al. 2006), but quasigeostrophic potential 
vorticity gradients may be more strongly negative on 
Saturn than they are in our simulation (Read et al. 2009). 
The vorticity profiles indicate that barotropic instability 
limits the sharpening of the retrograde jets, though not to 
the degree that the statistically steady states of the flows 
would satisfy sufficient conditions for linear barotropic 
stability for unforced and non-dissipative flows. (There 
is no reason that such sufficient conditions ought to be 
satisfied in forced-dissipative flows.) 

The jets are not only evident in long-term averages 
but also in instantaneous flow fields. The instantaneous 
zonal velocity and vorticity fields show the jets as well 
as large-scale jet undulations, waves, and coherent vor- 
tices (Fig. 3). In the equatorial region in the Jupiter and 
Saturn simulations, the waves are organized into large 
wave packets (Fig. 3, left column). Animations of the 
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Figure 2 Mean absolute vorticity (left) and quasigeostrophic 
potential vorticity (right) in the simulations, evaluated at the 
same levels at which the zonal velocities are shown in Fig. 1. 
The quasigeostrophic potential vorticity is calculated analo- 
gously to Read et al. (2009). It is not shown for the Jupiter 
simulation because the stratification at the corresponding level 
is nearly statically neutral (cf. Fig. 5 below), so potential vor- 
ticity is not well defined. 

flow fields^ show that the wave packets exhibit west- 
ward group propagation, as expected for long equatorial 
Rossby waves (Matsuno 1966; GiU 1982, chapter 11). In 
the vorticity fields (Fig. 3, right column), coherent vor- 
tices are seen in latitude regions with nearly homoge- 
nized absolute vorticity or potential vorticity, that is, in 
regions with large negative curvature of the zonal flow 
with latitude (cf. Fig. 2). 

b. High-latitude coherent vortices and waves on jets 

In high latitudes in the Jupiter and Saturn sim- 
ulations, very large coherent vortices (20° x 10° 
longitude X latitude) form spontaneously (Fig. 4). They 
extend all the way to the bottom of the domain, with the 
magnitude of the vorticity decreasing weakly with depth: 
the peak vorticity at the bottom of the domain is about 
80% of its maximum value in the column. 

The large coherent vortices are cyclonic, with typi- 
cal vorticities of magnitude ~2 x 10~^ s~^. They are 
advected by the flow in their environment and have a 
local temperature minimum at the center (~10 K lower 
temperature than the environment). These cyclonic vor- 
tices are long-lived, with life-spans apparently deter- 
mined by the radiative timescale (the timescale on which 
eddies can modify the mean flow). In the Jupiter simu- 
lation with an atmosphere of 3 bar thickness, the radia- 
tive timescale is ~10 Earth years; it is ~50 Earth years 

^available at http://www.gps.caltech.edu/~tapio/pubs.html 
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-160 160 ms"^ -4.0 0.0 4.0x10"V 




-240 240 ms-^ -3.0 0.0 3.0x10"V 




-50 50 ms-^ -2.0 0.0 2.0x10"^s"^ 




-50 50 ms-^ -2.0 0.0 2.0x1 0"^s"^ 

Figure 3 Zonal velocity (left column) and relative vorticity 
(right column) at one instant in the statistically steady state 
of the simulations. The levels at which the flow fields are 
shown are the same as in Fig. 1. The equatorial Rossby 
waves (organized into large wave packets) that are respon- 
sible for the generation of the equatorial superrotation are 
recognizable in the Jupiter and Saturn simulations. Coher- 
ent vortices are clearly seen in the Jupiter and Uranus sim- 
ulations. Animations of the flow fields are available at 
http://www.gps.caltech.edurtapio/pubs.htmL 

in the Saturn simulation. Since the radiative timescale 
increases with pressure, it is longer for deeper atmo- 
spheres, which might explain why the observed coher- 
ent vortices such as the Great Red Spot on Jupiter are so 
long-lived. 

Coherent vortices preferentially exist in regions where 



absolute vorticity or potential vorticity gradients van- 
ish, as they can then arise spontaneously in barotropic 
or quasigeostrophic flows and remain stable (e.g., 
McWilliams 1984; Marcus 1988, 1993). Since the plan- 
etary vorticity gradient vanishes at the poles, formation 
of coherent vortices in high latitudes may require less 
vorticity mixing in the environment than it does at lower 
latitudes. Hence, the large coherent vortices in high lat- 
itudes may appear earlier in simulations. If the simu- 
lations were conducted for a (much) longer period and 
if numerical (subgrid-scale) dissipation could be further 
reduced, it is possible that large coherent vortices would 
also appear in lower latitudes, such as the latitude (23° S 
planetocentric) of the Great Red Spot on Jupiter, which 
is embedded in an environment of small absolute vortic- 
ity and potential vorticity gradients (Ingersoll et al. 1981; 
Read et al. 2006). 

The coherent polar vortices in our simulations are con- 
tained in the polar cap bounded by the highest-latitude 
prograde jet. These polar jets exhibit large-scale undula- 
tions, as do other jets (cf. Fig. 3). In the Saturn simula- 
tion, for example, the prograde polar jet at 68° N exhibits 
a wavenumber-8 or 9 undulation (Fig. 4), with a zonal 
phase velocity of — 24ms~^. This is retrograde, and 
retrograde relative to the mean flow, consistent with the 
undulation being a Rossby wave. The undulation is rem- 
iniscent of the nearly stationary wavenumber-6 pattern 
("polar hexagon") observed in Saturn's polar atmosphere 
at 76°N planetocentric latitude (Godfrey 1988; Allison 
et al. 1990; Fletcher et al. 2008). Indeed, with a smaller 
drag coefficient that leads to a slightly stronger polar jet 
(see section 4d), we also obtain a wavenumber-6 pattern 
on the polar jet in our simulations; we will describe this 
in greater detail elsewhere. 

c. Vertical structure of zonal flow 

The simulated flows in Figs. 1-4 were shown near the 
suspected levels of observed cloud features on the giant 
planets. However, the flows in the simulations vary in 
the vertical. The prograde equatorial jets in the Jupiter 
and Saturn simulations strengthen with depth (Fig. 5, left 
column). The corresponding vertical shear of the zonal 
flow (~l-2 X 10~^s~^)is similar to that measured by 
the Galileo probe on Jupiter between 0.7 bar and 4 bar 
(Atkinson et al. 1998) and to that inferred from Cassini 
data for Saturn between 0.05 and 0.8 bar (Sanchez- 
Lavega et al. 2007; see also the zonal flow observations 
at different levels in Fig. 1). The retrograde equatorial 
jets in the Uranus and Neptune simulations are strongest 
in the stratosphere and weaken with depth, consistent 
with inferences drawn from gravity measurements with 
the Voyager 2 spacecraft (Hubbard et al. 1991). Away 
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Figure 4 Relative vorticity (10~^ s~^) in high latitudes at one 
instant in the statistically steady state of the Jupiter and Sat- 
urn simulations. Left column: south polar projection; right 
column: north polar projection. The vorticity is shown at the 
same levels as the flow fields in Figs. 1 and 3 (0.75 and 0.1 bar, 
respectively). 

from the equator, prograde jets generally weaken with 
depth and retrograde jets strengthen slightly or do not 
vary much with depth. 



d. Temperature structure 

Consistent with thermal wind balance, temperatures in- 
crease equatorward along isobars where prograde jets 
weaken with depth or retrograde jets strengthen with 
depth, and they decrease equatorward where the oppo- 
site is true. Therefore, in the equatorial upper tropo- 
sphere in the Jupiter and Saturn simulations, where the 
prograde jets strengthen with depth, temperatures de- 
crease equatorward and have a minimum at the equa- 
tor (Fig. 5, contours in right column). A similar equa- 
torial temperature minimum is seen in observations of 
Jupiter and Saturn (Simon-Miller et al. 2006; Fletcher 
et al. 2007). The tropopause, recognizable as the level at 
which the vertical temperature lapse rate changes sign, 
in all simulations lies near 0.1 bar, likewise as observed 
(Simon-Miller et al. 2006; Fletcher et al. 2007). Be- 
low the tropopause, temperatures increase with depth. In 
the Jupiter, Saturn, and Neptune simulations, the atmo- 
sphere is close to statically neutrally stratified below the 



statically stable layer near the tropopause because of vig- 
orous convection driven by intrinsic heat fluxes. In the 
Uranus simulation, the entire atmosphere is stably strati- 
fied (or close to it) because convection and intrinsic heat 
fluxes are weak (Fig. 5, colors in right column). 

In the Jupiter simulation, the equator-to-pole contrast 
in the brightness temperature of thermal radiation is 
10 K, which is similar to albeit larger than the observed 
brightness temperature contrast or the observed temper- 
ature contrast near the emission level (Ingersoll et al. 
1976; Ingersofl 1990; Ingersofl et al. 2004; Simon-Miller 
et al. 2006).^ In the Saturn simulation, the equator-to- 
pole brightness temperature contrast is 7 K, consistent 
with observations (Ingersoll 1990; Fletcher et al. 2007). 
In the Uranus and Neptune simulations, the equator- 
to-pole brightness temperature contrasts are 8 K and 
2 K, respectively, consistent with observations (Inger- 
soll 1990). That is, meridional enthalpy transport in 
all simulations substantially reduces the (much greater) 
radiative-convective equilibrium temperature contrasts 
near the emission levels. (The emission levels lie be- 
tween about 0.3 and 0.5 bar in our simulations, and 
radiative-convective equilibrium temperature contrasts 
there vary between 20 K for Neptune and 35 K for 
Uranus.) It does not appear necessary to invoke merid- 
ional mixing deep in the atmosphere to account for the 
smallness of the observed brightness temperature con- 
trasts (cf. Ingersoll 1976; Ingersoll and Porco 1978). 

However, although temperature contrasts at the emis- 
sion level are generally small, equator-to-pole temper- 
ature contrasts at higher or lower levels of the simu- 
lated atmospheres differ, and the same is likely true 
for the actual planets. For example, while temperature 
contrasts near the emission level in the Jupiter simula- 
tion are small, they are greater at lower levels where 
temperatures are greater (Fig. 6a). Because the atmo- 
sphere is close to statically neutrally stratified below 
the upper troposphere, entropy (potential temperature) 
there is constant in the vertical. [More generally, en- 
tropy is constant along angular momentum surfaces, to 
achieve a state of neutrality with respect to slantwise 
convection; see Emanuel (1983) and Thorpe and Ro- 
tunno (1989).] Hence, the meridional potential tempera- 
ture distribution at lower levels is the same as that near 
the top of the neutrally stratified layer, except near the 
equator where the atmosphere has a weak positive static 
stability (Fig. 6b). But this implies that meridional tem- 
perature gradients off the equator increase with pressure, 
as temperature T and potential temperature 6 are related 



^Variations in brightness temperature gradients off the equator may 
be weaker in the Jupiter and Saturn simulations than they are in obser- 
vations, at least at some wavelengths (Ingersoll 1990; Ingersoll et al. 
2004). 




Figure 5 Mean flow fields in the latitude-pressure plane in the simulations. The left column shows zonal-flow profiles: Gray 
contours for zonal flow speeds between 5 and 30 m s~^, with a contour interval of 5 m ; black contours for zonal flow speeds 
of 35 m or above, with a contour interval of of 35 m s~^. Solid contours and red tones for prograde flow and dashed contours 
and blue tones for retrograde flow. The right column shows temperature (contours, contour interval 10 K) and buoyancy frequency 
N (colors). The thick green parts of the latitude axes in the left column mark the latitudes with nonzero drag. The thin green lines 
indicate the levels at which flow fields are shown in Figs. 1-4. 



by T = 0{p/pqY, where p is pressure and po a con- 
stant reference pressure. In particular, the signatures of 
the prograde off-equatorial jets weakening with depth 
(enhanced meridional temperature gradients) and of the 
retrograde off-equatorial jets strengthening slightly or 
not varying with depth (reduced or vanishing meridional 
temperature gradients) are visible at all levels (Figs. 6a 
and b). There are no observations of entropy or temper- 
ature distributions below the upper troposphere for the 
giant planets, but we expect them to behave similarly as 
in our simulations, for the reasons discussed in section 5c 
below. 

4. Mechanisms of jet formation 

Why are the flow and temperature structures in the gi- 
ant planet simulations so different? The fundamental 
reason lies in the different strengths of the differential 
radiative heating and the intrinsic heat flux, and in the 
different ways in which the two can lead to the gen- 



eration of the eddies that maintain the jets. Eddies in 
rapidly rotating atmospheres generally transport angu- 
lar momentum from their dissipation (breaking) region 
into their generation region (Held 1975; Andrews and 
Mclntyre 1976, 1978; Rhines 1994; Held 2000; Vallis 
2006, chapter 12). If they are preferentially generated in 
prograde jets, they lead to angular momentum transport 
from retrograde into prograde jets, which can maintain 
the jets against dissipation (e.g., Vallis 2006; O'Gorman 
and Schneider 2008). Such angular momentum transport 
from retrograde into prograde jets has indeed been ob- 
served on Jupiter and Saturn (Ingersoll et al. 1981; Salyk 
et al. 2006; Del Genio et al. 2007). A central question is, 
then, what kind of eddies can give rise to the angular mo- 
mentum transport required to spin up and maintain the 
jets? We have addressed this question more formally and 
in greater detail in SL09. Here we summarize some re- 
sults from that earlier paper and expand on some that are 
important for understanding the simulations presented in 
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Figure 6 Temperature (a) and potential temperature (b) at the 
0.75-bar (solid), 1.0-bar (dashed), and 2.5-bar (dash-dotted) 
levels in Jupiter simulation. The potential temperature = 
T{po /pY is evaluated with the reference pressure po = 1 bar. 

this paper. 

a. Off- equatorial jets 

Away from the equator, the differential radiative heat- 
ing of the upper atmospheres produces meridional tem- 
perature gradients, which are baroclinically unstable and 
lead to eddy generation. Eddy generation preferentially 
occurs in the troposphere in the baroclinically more un- 
stable prograde jets with enhanced temperature gradi- 
ents and enhanced prograde vertical shear (Fig. 5). It 
results in angular momentum transport from retrograde 
into prograde tropospheric jets (Fig. 7). This angular 
momentum transport maintains the off-equatorial tropo- 
spheric jets in the Jupiter and Saturn simulations against 
dissipation at depth. It has a baroclinic structure and, 
in the Jupiter and Saturn simulations, is in structure 
and magnitude consistent with observations (cf. Inger- 
soll et al. 1981; Salyk et al. 2006; Del Genio et al. 2007). 
The conversion rate of eddy to mean-flow kinetic en- 
ergy is of order 10~^ W in the upper tropospheres 
in the Jupiter and Saturn simulations — as observed. In 
the global mean, the conversion rates are 0.09 Wm~^ 
and 0.026 W in the Jupiter and Saturn simulations, 
respectively, implying that the conversion rate in either 
simulation is about 0.6% of the total energy uptake of the 



atmosphere. Consistent with a baroclinic eddy genera- 
tion mechanism, off-equatorial eddy angular momentum 
fluxes and jets disappear in a Jupiter control simulation 
in which baroclinic instability is suppressed by imposing 
insolation uniformly at the top of the atmosphere (SL09). 

The off-equatorial jets in the Uranus and Neptune sim- 
ulations are situated in the stratosphere and are broader 
than those in the Jupiter and Saturn simulations. They 
appear to be broader because only the longest waves 
generated at lower levels are able to reach the strato- 
sphere (Charney and Drazin 1961). Indeed, the eddy 
angular momentum flux divergence at the stratospheric 
jet cores is much weaker than that at the tropospheric jet 
cores in the Jupiter and Saturn simulations (Fig. 7). It 
is not even of a consistent sign at all times but exhibits 
considerable low-frequency variability, as indicated by 
the differences between the statistically identical hemi- 
spheres in the long-term (1500-day) averages shown in 
Fig. 7. The jets become very weak below the tropopause 
and, particularly in the Uranus simulation, give way to 
a tropospheric zonal flow with smaller meridional scales 
(Fig. 5). The structure of the stratospheric jets implies 
that they interact only weakly with the drag at the lower 
boundary, so only weak eddy angular momentum flux di- 
vergence is necessary to maintain them. The jets are pri- 
marily a manifestation of the thermal structure and of the 
vertical shear of the zonal flow implied by it. The latter 
do not exhibit smaller-scale variations because smaller- 
scale eddy transports of angular momentum and heat in 
the stratosphere are weak. 



b. Equatorial superrotation 



Near the equator, convection can penetrate into the up- 
per troposphere and can generate Rossby waves if the in- 
trinsic heat flux is strong enough to overcome the static 
stabilization of the atmosphere by the radiative heat- 
ing from above. Fluctuations in convective heating are 
primarily balanced by vertical motion and, at the level 
of the convective outflows, by horizontal divergence of 
mass fluxes, as in the tropics of Earth's atmosphere. That 
is, the dominant balance in the thermodynamic equation 
is the weak-temperature gradient (WTG) balance (Sobel 




Figure 7 Mean zonal velocities, mass flux streamfunction (contours) and divergence div(i^^^;^ a cos 0) of meridional eddy angular 
momentum fluxes (colors) in the simulations. The left column shows zonal-flow profiles (contours as in Fig. 5) and eddy angular 
momentum flux divergence (colors). The right column shows the mass flux streamfunction (contours) and the same eddy angular 
momentum flux divergence as in the left column (colors). The contouring for the streamfunction is logarithmic: black contours 
from ±1 to ±64 x lO^kgs"^; gray contours for absolute values greater than or equal to 128 x lO^kgs"^, with factors of 
2 separating contour levels. (Black contours that would be within gray contours are not shown.) Solid contours for positive 
streamfunction values (counterclockwise rotation) and dashed contours for negative streamfunction values (clockwise rotation). 
Some streamfunction contours are truncated at the bottom of the plotting domain (3 bar) because they close at higher pressures 
(the pressure at the bottom of the GCM domain can locally exceed 3 bar). The contouring for the eddy angular momentum flux 
divergence is likewise logarithmic, with the scale shown in the colorbar. As in Fig. 5, the thick green parts of the latitude axes mark 
the latitudes with nonzero drag. 



et al. 2001),^ 

Here, denotes the divergent horizontal flow com- 
ponent, Q = DO/Dt the diabatic heating rate, and 

^The WTG approximation holds where the Rossby number sat- 
isfies Ro = U/\fL\ > 1 and the Froude number satisfies Fr = 
UV/ (gH) <C 1, where C/ is a zonal velocity scale, V a meridional 
or eddy velocity scale, L a length scale of flow variations, and H the 
scale height (Charney 1963; SL09). In SL09, we showed that the WTG 
approximation holds within ~4° of the equator in Jupiter's upper tro- 
posphere. Analogous scale analysis suggests the WTG approximation 
holds within ~7° of the equator in Saturn's upper troposphere. For 
Uranus' and Neptune's upper tropospheres, no flow data are avail- 
able, but with the tropospheric velocity scales from our simulations, 
the WTG approximation holds within ~5° of the equator. 



S = —dpO the static stability; the subscript h on the dif- 
ferential operator V/^ signifies horizontal derivative. As 
discussed in Sardeshmukh and Hoskins (1988), fluctua- 
tions in the horizontal divergence are a source of equa- 
torial Rossby waves, and the fluctuating vorticity source 

R' = R-R, (2a) 

with 

R=-U^h-^x)-{''x-^h)Ca. (2b) 

can be taken to be the Rossby wave source. (The over- 
bar denotes the isobaric zonal and temporal mean and 
primes deviations therefrom.) Thus, fluctuations in con- 
vective heating lead to horizontal divergence fluctua- 
tions (1), which can generate equatorial Rossby waves 
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through stretching of absolute vorticity (a or advection 
of absolute vorticity by the divergent flow (2). Be- 
cause the planetary vorticity vanishes at the equator, the 
Rossby wave source R' typically has largest amplitude 
just off the equator, but it does not necessarily vanish at 
the equator or where the absolute vorticity vanishes be- 
cause absolute vorticity advection by the divergent flow 
may not vanish. 

Equatorial Rossby waves, organized into large-scale 
wave packets, are recognizable in the Jupiter and Sat- 
urn simulations in Fig. 3. In the Jupiter simulation, 
the energy-containing zonal wavenumber is ^^10, cor- 
responding to the wavenumber of the wave packet enve- 
lope. Waves with similar scales have also been observed 
on Jupiter (Allison et al. 1990). The retrograde tilt of the 
waves' phase lines away from the equator, clearly seen in 
the Jupiter simulation (Fig. 3), indicates that they trans- 
port angular momentum toward the equator (cf. Peixoto 
and Oort 1992, chapter 11). This is generally to be ex- 
pected for such convectively generated Rossby waves: 
they transport angular momentum toward the equatorial 
region because this is where they are preferentially gen- 
erated. The Rossby wave source R' owing to horizon- 
tally divergent flow has largest amplitude in the equa- 
torial region because only there will convective heat- 
ing fluctuations necessarily lead to horizontal divergence 
fluctuations on large scales; away from the equator, the 
WTG approximation (1) of the thermodynamic equation 
does not hold. [See SL09 (their Fig. 5) for a demonstra- 
tion that R' has largest amplitude near the equator.] The 
angular momentum transport toward the equatorial re- 
gion by convectively generated Rossby waves leads to 
equatorial superrotation if it is sufficiently strong and 
drag on the zonal flow is sufficiently weak (SL09). 

Convective Rossby wave generation near the equator 
is the key process responsible for superrotation in the 
Jupiter and Saturn simulations. In SL09, we demon- 
strated that without intrinsic heat fluxes and the con- 
vection they induce, a Jupiter simulation similar to the 
one here exhibits equatorial subrotation; the same is true 
for the Jupiter and Saturn simulations here. Therefore, 
we suggest convective Rossby wave generation is what 
causes the superrotating equatorial jets on Jupiter and 
Saturn. 

When convective Rossby wave generation produces 
equatorial superrotation, it produces a jet whose half- 
width Lg is similar to the scale of equatorial Rossby 
waves: the equatorial Rossby radius = ^Jcj with 
gravity wave speed c and planetary vorticity gradient ^. 
Vorticity mixing arguments give an estimate for the max- 
imum strength of the equatorial jet (Rhines 1994, SL09): 
If the end state of vorticity mixing is a state in which the 
absolute vorticity is homogenized across the equatorial 



jet in each hemisphere separately, with a barotropically 
stable jump at the equator, and if the jet half- width is 
similar to the equatorial Rossby radius L^, the jet speed 
at the equator will be 



To the extent that this bound is attained (at the level of 
maximum equatorial jet speed), the jet speed increases 
quadratically with the jet width. This is roughly consis- 
tent with observations of Jupiter and Saturn (though the 
maximum equatorial jet speed on Saturn is not known 
for lack of observations deeper in the atmosphere). It is 
also roughly consistent with our simulations, although 
a state of homogenized absolute vorticity in the equa- 
torial region of each hemisphere is not attained in the 
simulations. That is, the equatorial jet on Saturn may 
be stronger and wider than that on Jupiter because the 
gravity wave speed is larger. 

The flow configurations in the Jupiter and Saturn sim- 
ulations differ qualitatively from those in the Uranus and 
Neptune simulations because the relative strengths of 
baroclinic eddy generation away from the equator and 
convective Rossby wave generation near the equator dif- 
fer. In the Uranus simulation, the intrinsic heat flux is 
negligible, the atmosphere is stably stratified, and there 
is no substantial convective Rossby wave source near the 
equator. Consequently, the equatorial eddy kinetic en- 
ergy is weak, and the equatorial flow is retrograde. 

In the Neptune simulation, the intrinsic heat flux is 
strong enough that convection penetrates into the up- 
per troposphere. As in the other simulations, eddies can 
be generated by (a) baroclinic instability off the equa- 
tor induced by differential solar heating, or (b) convec- 
tive Rossby wave generation near the equator induced 
by the intrinsic heat flux. Eddies produced by these two 
mechanisms compete with each other in their contribu- 
tion to the angular momentum transport to or from low 
latitudes. Off-equatorial baroclinic eddy generation im- 
plies angular momentum flux convergence in the off- 
equatorial generation regions and divergence in lower 
latitudes, and hence a tendency toward retrograde equa- 
torial flow. Convective Rossby wave generation near the 
equator can lead to prograde equatorial flow, but in the 
Neptune simulation, the rms Rossby wave source R! in 
the equatorial region is much smaller than that in the 
Jupiter and Saturn simulations: the rms Rossby wave 
source R! in the upper troposphere near the equator is 
^10~^^ for Neptune but ^10"^° for Jupiter and 
Saturn. Convective Rossby wave generation and the as- 
sociated angular momentum flux convergence near the 
equator appear to be too weak to overcome the angu- 
lar momentum flux divergence in low latitudes that is 
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caused by eddies generated baroclinically away from the 
equator. As a consequence, the equatorial flow is retro- 
grade. This is corroborated by control simulations. 

c. Neptune control simulations 

We investigated the relative roles of baroclinic eddy gen- 
eration and convective Rossby wave generation on Nep- 
tune in two control simulations, one in which the convec- 
tive Rossby wave source was enhanced and one in which 
the baroclinic eddy generation caused by differential so- 
lar heating was suppressed. Because angular momen- 
tum flux divergence in low latitudes owing to baroclinic 
eddy generation away from the equator can counteract 
any angular momentum flux convergence owing to con- 
vective Rossby wave generation, generation of equato- 
rial superrotation in the Neptune simulation may require 
a stronger intrinsic heat flux or weaker differential solar 
heating. 

Indeed, a control simulation with Neptune's physical 
parameters but in which the convective Rossby wave 
source was enhanced by enhancing the intrinsic heat 
flux — setting it to Saturn's 2.01 W in place of Nep- 
tune's 0.433 Wm~^ — exhibits equatorial superrotation 
(Fig. 8, left column). Conversely, a control simulation 
in which baroclinic eddy generation was suppressed by 
imposing insolation uniformly at the top of the atmo- 
sphere (but keeping the global mean fixed) also exhibits 
equatorial superrotation (Fig. 8, right column). The pro- 
grade off-equatorial jets disappear with the suppression 
of baroclinicity. 

We suggest, then, that Uranus and Neptune exhibit 
equatorial subrotation because baroclinic eddy genera- 
tion away from the equator is strong compared with con- 
vective Rossby wave generation near the equator. Inter- 
estingly, our simulations suggest that Neptune's intrinsic 
heat flux only needs to be larger by an 0(1) factor for 
Neptune's atmosphere to develop equatorial superrota- 
tion, implying that Neptune may have been superrotat- 
ing earlier in its history, when intrinsic heat fluxes were 
stronger. 

d. Effect of drag formulation on simulated flows 

One relatively unconstrained aspect of our simulations is 
the strength and functional form of the drag at the arti- 
ficial lower boundary. We investigated the sensitivity of 
our results to the drag formulation by varying it in a few 
Jupiter simulations. 

The Jupiter simulation in SL09 had an equatorial no- 
drag region half as wide as that here (extending to (po = 
16° latitude vs. 0o = 33° here), in addition to having a 
larger drag coefficient off the equator [ko = 1/(20 days) 



vs. ko = 1/(100 days) here]. Mean flow fields from that 
earlier simulation are shown in Fig. 9a. The half- width 
of the prograde equatorial jet is ^5° smaller than in the 
Jupiter simulation reported here (cf. Fig. 5). The adja- 
cent strong retrograde jets appear to be confined to the 
no-drag region and hence do not extend as far poleward 
as in the simulation with the wider no-drag region. The 
off-equatorial jets are somewhat weaker and narrower — 
a result of the enhanced off-equatorial drag, consistent 
with theories and other simulations of geophysical tur- 
bulence (Smith et al. 2002; Danilov and Gurarie 2002). 
However, the width of the equatorial no-drag region does 
not primarily control the strength or width of the equato- 
rial jet, as evidenced by the relatively moderate changes 
in the flow in low latitudes in response to the factor 2 
change in the width of the no-drag region. 

That the strength and width of off-equatorial jets de- 
pend on the drag coefficient is directly illustrated by sim- 
ulations in which we increased the off-equatorial drag 
coefficient further [Fig. 9b, ko = 1/(10 days)]. The 
off-equatorial jets become weaker and narrower as the 
drag coefficient is increased. However, if the same en- 
hanced drag is used at all latitudes, without an equato- 
rial no-drag region, there is no large-scale prograde jet 
at the equator, while the off-equatorial flow is not sub- 
stantially modified (Fig. 9c). (But a narrow and shallow 
prograde jet forms near the tropopause at the equator.) 
Similarly, if a weaker constant drag is used at all lati- 
tudes [ko = 1/(100 days)], without an equatorial no- 
drag region, there is likewise no large-scale prograde jet 
at the equator (Fig. 9d). Even weaker equatorial drag is 
required to obtain a large-scale prograde jet at the equa- 
tor if intrinsic heat fluxes are specified consistent with 
observations. Consistent with the theoretical arguments 
in SL09, stronger equatorial drag requires larger intrinsic 
heat fluxes and thus a stronger equatorial Rossby wave 
source to lead to superrotation. However, neither the pre- 
cise functional form of the drag, nor the magnitude of the 
drag coefficient where it is nonzero, nor the width of the 
no-drag region appear to be essential for our results — 
as long as there is an equatorial region with no or suffi- 
ciently low drag such that a large-scale prograde jet can 
form.^ 

The simulations with different drag formulations show 
that better fits to observations can be obtained if different 
drag formulations are used for the different giant planets. 
This is physically justifiable because the interior proper- 
ties of the planets differ and give rise to differences in 
the strength of MHD drag and in the depth at which it 
acts (Liu et al. 2008). 

^The same arguments may also explain the formation of prograde 
equatorial jets in Scott and Polvani's (2008) simulations with an essen- 
tially frictionless shallow- water model. 
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Figure 8 Mean zonal velocities in the latitude-pressure plane in Neptune control simulations. Contour intervals and colors as in 
Fig. 5. Left: simulation with Neptune's physical parameters but Saturn's intrinsic heat flux (2.01 W m~^). Right: simulation with 
Neptune's physical parameters but uniform insolation at the top of the atmosphere. 




Figure 9 Mean flow fields in the latitude-pressure plane in Jupiter simulations with different drag formulations. The left column 
shows the zonal flow and the right column the temperature and buoyancy frequency, with the same plotting conventions and contour 
intervals as in Fig. 5. (a) No drag in the equatorial region up to = 16°, and off-equatorial drag coefficient ko = 1/(20 days), 
(b) No drag in the equatorial region up to = 33°, and off-equatorial drag coefficient ko = 1/(10 days), (c) Constant drag at all 
latitudes with drag coefficient ko = 1/(10 days), (d) Constant drag at all latitudes with drag coefficient ko — 1/(100 days). 



5. Mean meridional circulations and angular mo- 
mentum balance in deep atmospheres 

Our theory and simulations are consistent with the en- 
ergy and angular momentum balances of the giant plan- 
ets as far as they are known, and they are broadly 



consistent with many observed upper-atmospheric flow 
features. Their relevance, however, depends on how 
the flows in the upper atmospheres couple to flows at 
depth. We have represented this coupling in an ideal- 
ized fashion in our thin-shell simulations through the 
drag formulation. Here we show how results for the 
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upper-tropospheric flows constrain the flows at depth. 
What follows is a straightforward generalization of well 
known results for thin atmospheres — particularly the 
principle of "downward control" (Haynes et al. 1991) — 
which was already sketched in SL09. We give the argu- 
ments in some detail, as their implications for planetary 
atmospheres are underappreciated. 

a. Local angular momentum balance 

In any atmosphere, regardless of its constitutional law, 
the balance of angular momentum around the planet's 
spin axis can be written as 

dt{pM) + V • (puM) = -dxv + r^pV, (4) 

where M = M^q + Mu is the angular momentum per 
unit mass, composed of the planetary angular momen- 
tum Mq^ = Q.r\_ and the relative angular momentum 
Mu = ur±_ . Here, r_L = r cos (j) is the (cylindrically 
radial) distance to the planet's spin axis and r the (ra- 
dial) distance to the planet's center; P is a zonal drag 
force per unit mass, which may include viscous dissipa- 
tion; A is longitude (azimuth); and p is pressure and p 
density (e.g., Peixoto and Oort 1992, chapter 11). In the 
thin-shell approximation, the distance to the spin axis is 
approximated as = a cos 0, with a constant planetary 
radius a. But the angular momentum balance (4) also 
holds in a deep atmosphere if r_L = r cos is taken to be 
the actual distance to the spin axis, with variable r. 

In a statistically steady state, upon averaging tempo- 
rally and zonally (azimuthally), the angular momentum 
balance becomes 

u* • VMn + u* • VM* = vaV" - S, (5) 

where 

S=-_V- (pS^*) (6) 

is the eddy angular momentum flux divergence. The 
overbar (•) now denotes the temporal and zonal mean at 

constant r_L, and (•) = {p •) /p denotes the correspond- 

— * 

ing density- weighted mean; primes {-)' = (•) — (•) 
denote deviations from the latter. The eddy angular 
momentum flux divergence in Fig. 7 is the pressure- 
coordinate analog of the meridional component of the 
flux divergence (6). 

The ratio of the second to the first term on the left- 
hand side of the angular momentum balance (5) is of 
order Rossby number. 



where is the length scale of flow variations in the 
cylindrically radial direction. That is, if L is a merid- 
ional length scale, L±_ = L sin ^ is the projection of the 
meridional length scale onto the equatorial plane, and 
Ro = U /\fL\ becomes the familiar Rossby number for 
the thin- shell approximation. Away from the equator, the 
Rossby number is generally small in the tropospheres 
of the giant planets if zonal flow velocities at depth do 
not substantially exceed those observed on Jupiter and 
Saturn, or those seen in the tropospheres of Uranus and 
Neptune in our simulations (see SL09 and footnote 5). 
The angular momentum balance then is approximately 

u* • VMn ^ r_Lp* - S. (8) 

The term on the left-hand side represents the advection 
of planetary angular momentum by the mean flow, or the 
Coriolis torque per unit mass (u* • VMq = — /^* r_L 
in the thin-shell approximation). Three special dominant 
balances can be distinguished. 

1). V^O,Sj^O 

This is the dominant balance in the off-equatorial upper 
troposphere, where eddy angular momentum flux diver- 
gences are significant but drag forces are negligible. In 
this case, the angular momentum balance 

u* . VMn « -S (9) 

implies that the mean mass flux has a component across 
Mq surfaces: toward the planet's spin axis (poleward) 
where eddy angular momentum fluxes diverge (S > 0), 
and away from the spin axis (equatorward) where they 
converge (S < 0). Because eddy angular momen- 
tum fluxes in our simulations generally diverge in ret- 
rograde tropospheric jets, or zonal velocity minima, and 
converge in prograde tropospheric jets, or zonal veloc- 
ity maxima, the mean meridional mass flux in the off- 
equatorial upper troposphere is generally poleward in 
retrograde jets and equatorward in prograde jets (Fig. 7). 
The same is almost certainly true on Jupiter and Sat- 
urn, where similar eddy angular momentum fluxes in 
the upper troposphere have been observed (Salyk et al. 
2006; Del Genio et al. 2007). (Eddy angular momentum 
fluxes have not been observed on the other giant planets.) 
Direct observations of mean meridional mass fluxes on 
Jupiter and Saturn are ambiguous, but the observed dis- 
tribution of convection provides indirect evidence that 
there is upwelling in the cyclonic shear zones between 
retrograde and prograde jets (Ingersoll et al. 2000; Porco 
et al. 2003; Del Genio et al. 2007), consistent with these 
arguments and with our simulations (Fig. 7). 
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2) . V^O,S^O 

This is the dominant off-equatorial balance immediately 
below the layer with significant eddy angular momentum 
flux divergences, where drag forces are negligible. In 
this case, the angular momentum balance 

u* • VMn ^ (10) 

implies that the mean mass flux is along Mq surfaces, 
that is, parallel to the planet's spin axis in deep atmo- 
spheres or vertical in thin atmospheres. As in Earth's at- 
mosphere, such an off-equatorial tropospheric layer with 
mean mass flux along Mq surfaces is clearly seen in 
our simulations, where Mq surfaces are vertical (Fig. 7). 
It very likely also exists at least on Jupiter and Saturn, 
where, as we argued in the introduction, energetic con- 
straints indicate that significant eddy angular momentum 
fluxes cannot extend deeply into the atmosphere. 

The constraint that the mean mass flux is along Mq 
surfaces is not to be confused with the Taylor-Proudman 
constraint. The Taylor-Proudman constraint states that 
steady-state velocities in rapidly rotating barotropic at- 
mospheres do not vary in the direction of the planet's 
spin axis if non-conservative forces are absent (e.g., 
Kaspi et al. 2009). It requires the flow to be barotropic, 
whereas the flows we consider generally are baroclinic 
and sheared along Mq surfaces, as in Earth's atmosphere 
and in our simulations (e.g.. Figs. 5 and 6). 

3) . Vj^O,S^O 

This is the dominant off-equatorial balance (Ekman bal- 
ance) in lower layers in our simulations, where drag 
forces are significant. It very likely also is the dominant 
balance in any deep layer of significant drag on the giant 
planets. In this case, the angular momentum balance 

u* -VMn^VAV" (11) 

implies that the mean mass flux has a component across 
Mq surfaces: toward the planet's spin axis (poleward) 
where the drag force is retrograde (V < 0), and away 
from the spin axis (equatorward) where it is prograde 
(V > 0). To the extent that the drag force locally retards 
the mean zonal flow, as it does for the linear drag in our 
simulations, it implies that away from the equator, there 
is a mean mass flux toward the spin axis where the mean 
zonal flow is prograde and away from the spin axis where 
it is retrograde (Fig. 7). 

b. Mean meridional circulation and zonal flow at depth 

Thus overturning mass circulations in the meridional 
plane come about. In a statistically steady state, any 



mean mass flux across an Mq surface associated with 
eddy angular momentum flux divergences in the upper 
troposphere must be balanced by an equal and opposite 
mean mass flux across the same Mq surface somewhere 
else, to obtain closed circulation cells. Where the Rossby 
number is small, this opposing mean mass flux must be 
associated with an opposing eddy angular momentum 
flux divergence or drag. Outside the equatorial no-drag 
region in our simulations, the opposing mean mass flux 
is associated with drag at depth, similar to how mass cir- 
culation cells close in Earth's atmosphere. On the giant 
planets, MHD drag acts at great depth and can fulfill a 
similar role in closing circulation cells. 

The angular momentum balance also constrains the 
zonal flow at depth. Taking a density- weighted integral 
of the angular momentum balance (8) along Mq surfaces 
and using mass conservation shows that any net diver- 
gence or convergence of eddy angular momentum fluxes 
on an Mq surface must be balanced by a zonal drag force 
on the same Mq surface, 

{pS}n^r^{pV*}n, (12) 

where {'}q denotes an average over Mq surfaces. To 
the extent that the drag force locally retards the mean 
zonal flow, it follows that if eddy angular momentum 
flux convergence occurs in the upper troposphere in pro- 
grade jets, and divergence in retrograde jets, and if this 
is the dominant eddy angular momentum flux conver- 
gence/divergence on an Mq surface, the mean zonal flow 
where the drag acts must be of the same sign as the flow 
in the upper troposphere on the same Mq surface. That 
is, zonal jets must extend to wherever drag acts, irre- 
spective of its depth, even if the eddy angular momen- 
tum fluxes are confined to the upper troposphere [see 
O' Gorman and Schneider (2008) for a numerical exam- 
ple]. Because drag cannot act at the upper boundary of 
the atmosphere (it would imply an impossible torque on 
outer space), the jets generally extend downward. In- 
sofar as the eddy angular momentum fluxes in the up- 
per atmosphere control the dissipation at depth, one may 
speak of "downward control" of the mean meridional cir- 
culation and zonal flow at depth (Haynes et al. 1991). 
Figure 10 summarizes these inferences from the off- 
equatorial angular momentum balance. 

The mean meridional circulation cells link the dynam- 
ics in the upper troposphere to the flow at depth, ad- 
justing the thermal structure of the atmosphere between 
the upper troposphere and the layer where the drag acts 
such that the zonal-flow shear along Mq surfaces im- 
plied by thermal- wind balance becomes consistent with 
the balance (12) between eddy angular momentum flux 
divergences and drag. This is analogous to the role mean 
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^ Prograde flow 



^ Retrograde flow 




Figure 10 Sketch of mean meridional circulation and zonal 
flow off the equator in giant planet atmospheres. Straight blue 
lines with arrows indicate the mass circulation; green lines in- 
dicate Mq contours; wavy lines indicate eddy angular momen- 
tum fluxes. The size of the zonal-flow symbols is to suggest the 
speed of the flow. The blue shaded region represents the elec- 
trically conducting part of the atmosphere, where MHD drag 
acts. 

meridional circulations play in Earth's atmosphere (e.g., 
Holton 2004, chapter 10). 

Where the Rossby number is not small, the advec- 
tion of angular momentum by the mean mass flux — the 
second term on the left-hand side of (5) — cannot be ne- 
glected and also contributes to the angular momentum 
balance and its density-weighted integral along Mq sur- 
faces, as discussed in SL09. However, in the tropo- 
spheres of the giant planets, this appears to be signif- 
icant only within a few degrees latitude of the equator 
(see SL09 and footnote 5). 

c. Implications for thermal structure 

Jupiter, Saturn, and Neptune have sufficiently strong in- 
trinsic heat fluxes to lead to convection in their tropo- 
spheres (Guillot et al. 2004; Guillot 2005), as in our sim- 
ulations. The occurrence of convection and the expected 
homogenization of entropy along convective plumes fur- 
ther constrains the thermal structure of the tropospheres 
and thus, through thermal wind balance, the zonal flow 
structures. 

In the giant planet tropospheres, the convective 



Rossby number is generally small, and viscous mo- 
mentum dissipation and thermal diffusion are thought 
to be negligible. Under these circumstances, convec- 
tive plumes are columns aligned with Mq^ surfaces, be- 
cause, as above, rapid rotation inhibits motion perpen- 
dicular to surfaces in the absence of viscous or 
other stresses (e.g., Busse 1976, 1978; Kaspi et al. 2009; 
Jones et al. 2009).^ Therefore, convection tends to ho- 
mogenize entropy along columns in the direction of the 
planet's spin axis (in deep atmospheres) or in the ver- 
tical (in thin atmospheres); however, it does not con- 
strain entropy gradients in perpendicular directions. The 
forced-dissipative statistically steady state that results in 
the presence of vigorous convection thus is neutral with 
respect to slantwise convective instabilities, that is, con- 
vective and inertial axisymmetric instabilities (Emanuel 
1983; Thorpe and Rotunno 1989; Emanuel 1994, chap- 
ter 12; Schneider 2007). 

In our Jupiter, Saturn, and Neptune simulations, such 
a state with nearly convectively neutral interior tropo- 
spheres is indeed attained outside a few degrees of the 
equator (Fig. 5). The entropy and its meridional gradient 
hence do not vary in the vertical (e.g.. Fig. 6), and neither 
does, by thermal wind balance, the shear (with respect to 
pressure) of the zonal flow. A corresponding state with 
entropy and its meridional gradient homogeneous in the 
direction of the spin axis can be expected to be attained 
on the actual planets, and thermal wind balance for deep 
atmospheres then similarly constrains the zonal flow. For 
a deep atmosphere in the anelastic approximation (valid 
for small fluctuations relative to an isentropic reference 
state), thermal wind balance reads 

oz r 0(p 

where z is the cylindrical height coordinate in the di- 
rection of the spin axis, s is the specific entropy, and 
as{p) is a thermal expansion coefficient that relates en- 
tropy fluctuations to density fluctuations (Kaspi et al. 
2009). That is, if the meridional entropy gradient does 
not vary in the direction of the spin axis, the shear of 
the zonal flow in the direction of the spin axis depends 
only on a^, g, and r, all of which generally vary (as 
and g vary primarily with r). In Jupiter's and Saturn's 
upper tropospheres, the meridional entropy gradient and 
thus the zonal thermal wind shear approximately van- 
ish at the zeros of the zonal wind (Simon-Miller et al. 
2006; Read et al. 2006, 2009), as they do in our sim- 
ulations (Fig. 5). To the extent that entropy deeper in 
the troposphere is homogenized in the direction of the 

^More generally, where the convective Rossby number is not nec- 
essarily small, convective plumes are aligned with angular momentum 
(M) surfaces. 
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spin axis, the thermal wind balance (13) suggests that 
the zonal flow shear then vanishes at all depths extend- 
ing downward in the direction of the spin axis from the 
upper-tropospheric zeros of the zonal flow. So zeros of 
the zonal flow project downward along cylinders con- 
centric with the spin axis.^ 

In the literature on the giant planets, it is often taken 
as axiomatic that their interior tropospheres are rendered 
isentropic by convection, resulting in zonal flows (Taylor 
columns) without shear in the direction of the spin axis 
(e.g., Vasavada and Showman 2005). However, convec- 
tion in general does not homogenize entropy in direc- 
tions perpendicular to angular momentum surfaces, so 
an isentropic interior cannot be assumed a priori. The 
flow on the giant planets, where the Rossby number is 
small, must satisfy the constraints (8)-(12) dictated by 
the angular momentum balance, as well as thermal wind 
balance (13). Given that significant eddy angular mo- 
mentum flux divergences in the upper tropospheres have 
been observed at least on Jupiter and Saturn, it is very 
unlikely that the zonal drag, which depends on the zonal 
flow, can balance the net eddy angular momentum flux 
divergence on an surface without any shear of the 
zonal flow in the direction of the spin axis. It hence is 
very unlikely that the interiors of the giant planets are 
isentropic.^ 

d. Implications for role of drag 

Significant eddy angular momentum flux divergences 
have been observed in Jupiter's and Saturn's upper tro- 
pospheres, but no mechanism has been proposed of how 
they could be balanced by opposing eddy angular mo- 
mentum flux divergences at depth such the angular mo- 
mentum balance (12) integrated over Mq^ surfaces is sat- 
isfied without a drag mechanism. On the other hand, 
coupling of the flow at depths at which the atmosphere 
is electrically conducting to the magnetic field is a plau- 
sible mechanism that generates MHD drag (Liu et al. 
2008; see also the appendix). This MHD drag can close 
the angular momentum balance integrated over sur- 
faces. Therefore, we adopted as our working hypothesis 
that the MHD drag acting at depth couples to the flow 
in the upper atmosphere, although details such as how 
the MHD drag depends on the zonal flow are poorly un- 

^ Where latent heat release in phase changes of water is dynamically 
important, a moist entropy rather than a dry entropy should be consid- 
ered, and such a moist entropy can be expected to be homogenized in 
the direction of the spin axis. 

^If entropy deviations from an isentropic reference state are not 
small so that the anelastic approximation cannot be made, the ther- 
mal wind balance (13) becomes more complicated (Kaspi et al. 2009). 
But this does not affect the quahtative considerations on which our 
conclusions are based. 



derstood. This approach gave statistically steady states 
in which the angular momentum balance is closed in a 
manner that is physically plausible and consistent with 
observations. 

Outside a few degrees latitude of the equator, the 
upper-tropospheric dynamics are linked to the flow and 
any drag at depth along cylindrical Mq^ surfaces. There- 
fore, outside the tangent cylinder that just grazes the re- 
gion of substantial MHD drag in the equatorial plane, 
the upper-tropospheric flow cannot be linked to drag at 
depth. We chose to represent this equatorial region of no 
effective drag on the upper-tropospheric flow in our thin- 
shell simulations by having an equatorial no-drag region. 
A no-drag region extending to 33° latitude corresponds 
to assuming that the region of substantial MHD drag is 
confined within about cos (33°) =0.84 planetary radii. 
It is doubtful that this is an accurate estimate for all giant 
planets (Liu et al. 2008). However, as we have shown in 
section 4d, where exactly the MHD drag acts, and how 
strong it is, does not affect the essence of our results. 

The preceding discussion shows that simulations of 
only thin atmospheric shells can have closed energy and 
angular momentum balances that are physically plausi- 
ble and consistent with observations of the giant plan- 
ets. A model domain of sufficient depth to take into ac- 
count the absorption of solar radiation in the upper at- 
mosphere is essential to obtain baroclinic flows with an 
energy balance that is consistent with observations. Re- 
solving baroclinic eddy fluxes of angular momentum in 
the upper atmosphere is essential to obtain an angular 
momentum balance that is consistent with observations. 
Drag at depth closes the balances in a physically plau- 
sible manner (e.g., without assuming excessive viscous 
stresses in the planetary interior). However, the depth of 
the nearly inviscid interior layer in which there are no 
significant heat sources and where P ^ and 5^0 
is immaterial for the mechanisms we discussed. The 
depth of that layer can be expected to affect quantitative 
aspects such as the vertical shear of the zonal flow, but 
we do not expect it to affect the qualitative aspects and 
large-scale flow features on which we have focused. 

6. Conclusions 

We have presented the first simulations of all four gi- 
ant planets with closed energy and angular momentum 
balances that are consistent with observations. The sim- 
ulations reproduce many large-scale features of the ob- 
served flows, such as equatorial superrotation on Jupiter 

^^Latent heat release in phase changes of water may play a role in 
that layer, but it does not represent an external heat source, merely a 
conversion between forms of energy, and hence it does not affect the 
integrated energy balance. 
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and Saturn and equatorial subrotation on Uranus and 
Neptune. They exhibit temperature structures that are 
broadly consistent with available observations, and they 
reproduce many details of the observed flows, for exam- 
ple, their vertical structure to the extent it is known and 
characteristic equatorial waves observed on Jupiter. We 
have demonstrated that equatorial superrotation is gener- 
ated if convective Rossby wave generation is strong and 
low-latitude angular momentum flux divergence owing 
to baroclinic eddies generated off the equator is suffi- 
ciently weak (Jupiter and Saturn); equatorial subrotation 
results if either convective Rossby wave generation is 
weak or absent (Uranus) or low-latitude angular momen- 
tum flux divergence owing to baroclinic eddies is suffi- 
ciently strong (Neptune). 

Current computational resources limit our ability to 
simulate flows at depth. However, considerations of the 
angular momentum balance have shown that the zonal 
jets should extend — generally with shear — to the depth 
where drag acts on them and balances the angular mo- 
mentum flux divergences and convergences in the upper 
troposphere. That drag acts on the zonal flow at depth is 
suggested by observations of eddy angular momentum 
fluxes on Jupiter and Saturn, and a plausible MHD drag 
mechanism exists. Though quantitative aspects (e.g., jet 
strength and shear) may be affected by our inability to re- 
solve the flow and drag at depth, the jet formation mech- 
anisms we discussed are not affected by it. 

We expect as-yet unobserved aspects of the flow and 
temperature structures to be consistent with the simula- 
tions and mechanisms we presented. For example, we 
predict that NASA's upcoming JUNO mission to Jupiter 
will find evidence of zonal flows with vertical shear simi- 
lar to those in Fig. 5: near the equator, a strong and deep 
prograde jet, and away from the equator, prograde jets 
that weaken and retrograde jets that weaken only slightly 
or strengthen with depth. The thermal and gravitational 
signature of such zonal flows will likely be measurable 
by JUNO. 
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APPENDIX A 

General Circulation Model 



a. Resolution 

The GCM solves the hydrostatic primitive equations us- 
ing the spectral transform method in the horizontal and 
finite differences in the vertical. The horizontal spectral 
resolution depends on the radius of the planet being sim- 
ulated (Table 1). The vertical coordinate is a = p/ps 
(pressure p normalized by pressure at lower boundary 
Ps); it is discretized with 30 levels for Jupiter and Saturn 
and 40 levels for Uranus and Neptune. 

b. Drag at lower boundary 

All parameter choices are constrained by knowledge of 
the physical properties of the planets and material prop- 
erties of their atmospheres, as well as by observations 
where available. However, the drag formulation at the ar- 
tificial lower boundary of the GCM is poorly constrained 
by data or physics. It represents the MHD drag the flow 
experiences in the interior of the planets. 

In the interior of Jupiter and Saturn, the conductiv- 
ity of molecular hydrogen increases with depth and be- 
comes approximately constant where hydrogen becomes 
metallic at ^1.4 Mbar (Nellis et al. 1996). In the interior 
of Uranus and Neptune, the conductivity of the gas enve- 
lope likewise increases with depth and is determined by 
the conductivity of hydrogen and water ice (Nellis et al. 
1997). In the high-conductivity interior, the interaction 
of the magnetic field and the fluid flow produces Ohmic 
dissipation and retards the flow (Liu et al. 2008). 

We represented this MHD drag deep in the atmo- 
sphere in the simplest possible way in our thin-shell 
GCM, choosing the same drag formulation and depth 
of the artificial lower boundary in all giant planet sim- 
ulations, to rule out that differences among them are 
caused by differences in poorly constrained parameters. 
We assume linear (Rayleigh) drag acts near the GCM's 
lower boundary, but only outside an equatorial latitude 
band (see SL09 and section 5d). As in the models in 
SL09 or Held and Suarez (1994), the drag coefficient 
decreases linearly from its value ko at the lower bound- 
ary at a = 1.0 to zero above a = 0.8. The equato- 
rial no-drag region extends to 00 = 33° latitude in all 
our simulations, corresponding to a MHD drag that acts 
only within 0.84 planetary radii. The drag coefficient is 
constant (ko = 10~^days~^) outside this region. The 
kinetic energy dissipated by the Rayleigh drag (a few 
percent of the sum of the intrinsic heat flux and the ab- 
sorbed solar radiative flux) is returned to the flow locally 
as heat to conserve energy. 

We chose the width of the no-drag region and the drag 
coefficient outside of it empirically, to obtain jets in the 
upper atmosphere that have similar strength and width 
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as the observed jets. By choosing drag formulations that 
differ from planet to planet, better fits to observations 
could be obtained (cf. section 4d). 

c. Radiative transfer 

As in SL09, radiative transfer is represented as that in 
a homogeneous gray atmosphere, using the two-stream 
approximation. The top-of-atmosphere (TOA) insola- 
tion is imposed as perpetual equinox with no diurnal cy- 
cle, Ftoa = (Fo/7r)cos0, with the appropriate solar 
constant Fq for each planet (Table 1). That is, for the 
purposes of this paper, we ignore the seasonal cycle in 
TOA insolation, which is substantial for Saturn, Uranus, 
and Neptune because of their obliquities. Ignoring sea- 
sonality may be justifiable, for example, if response 
timescales of the atmospheric circulation (e.g., radiative 
timescales) are much longer than seasonal timescales, 
which may be the case on the giant planets. However, 
the nonzero obliquities also influence the annual-mean 
TOA insolation, especially on Uranus, an influence we 
likewise ignore. 

The solar radiative flux for a semi-infinite scattering 
and absorbing atmosphere is calculated for a solar opti- 
cal depth Ts that is linear in pressure, Tq = Tso{p/po), 
to represent scattering and absorption by a well-mixed 
absorber. The solar optical properties of the giant planet 
atmospheres are not well constrained. To minimize dif- 
ferences among the simulations, we chose the same so- 
lar optical properties for all giant planets: Tso = 3.0 at 
Po = 3.0 bar. This gives a solar radiative flux quali- 
tatively consistent with Galileo probe measurements in 
Jupiter (Sromovsky et al. 1998). 

The thermal radiative flux is calculated for a gray 
atmosphere in which the thermal optical depth ti is 
quadratic in pressure, ri = r^o(p/Po)^. to represent 
collision-induced absorption of thermal radiation. The 
thermal optical depths t/q at pressure po are chosen such 
that the observed thermal emission levels (e.g., IngersoU 
1990) of the giant planets approximately correspond to 
n = 1. The thermal optical depths thus vary from planet 
to planet (Table 1). 

The lower boundary condition for the radiative fluxes 
is energy conservation: the upward thermal radiative flux 
is set equal to the sum of the downward solar flux and 
thermal radiative flux at each grid point. 

d. Intrinsic heat flux 

A spatially uniform and temporally constant heat flux, 
corresponding to that estimated for the giant planets (Ta- 
ble 1), is deposited in the GCM's lowest layer to mimic 
intrinsic heat fluxes. 



e. Convection scheme 

A quasi-equilibrium convection scheme represents (dry) 
convection. It relaxes temperature profiles toward a 
convective profile with adiabatic lapse rate F = g/cp 
(Schneider and Walker 2006). The convective relaxation 
time is chosen to be roughly the time it takes a gravity 
wave with speed c to traverse the extratropical Rossby 
radius = c//, that is, roughly an extratropical iner- 
tial time f~^. We chose the convective relaxation time to 
be 6 hr for Jupiter and Saturn and 10 hr for Uranus and 
Neptune. We experimented with convective relaxation 
times up to a factor 2 smaller and a factor 4 larger in pre- 
liminary simulations; the simulated flows appeared not 
to be sensitive to such variations of the relaxation time. 

The convection scheme does not transport momen- 
tum, that is, it assumes a convective Prandtl number 
of zero. Prandtl numbers for dry convection are usu- 
ally greater than zero, so this represents an idealization, 
which may affect our results. However, the convective 
generation of Rossby waves on which we focus occurs 
in the equatorial upper troposphere, where the vertical 
shear is relatively small (Fig. 5). So one may expect 
convective momentum fluxes not to alter our results sub- 
stantially. But to the extent that the vertical shear and 
convective momentum fluxes cannot be neglected, they 
may amplify the superrotation in the upper troposphere 
in the Jupiter and Saturn simulations, as the equatorial 
zonal flow in those simulations is stronger at depth than 
in the upper troposphere, so that convection can be ex- 
pected to transport momentum upward. The dependence 
of our results on the convective Prandtl number is worth 
investigating further. 

/ Subgrid-scale dissipation 

For a < 0.8, above the layer with Rayleigh drag, hor- 
izontal hyperdiffusion in the vorticity, divergence, and 
temperature equations is the only frictional process. The 
hyperdiffusion is represented by an exponential cutoff 
filter (Smith et al. 2002), with a damping time scale of 
2 h at the smallest resolved scale. The cutoff wavenum- 
ber depends on the horizontal resolution (Table 1). 

The energy dissipated by the subgrid-scale hyperdif- 
fusion is not returned to the flow as heat. However, it 
amounts to less than 1% of the total energy uptake of the 
atmosphere in all simulations. 

g. Simulations 

The simulations were spun-up from an isothermal rest 
state, with small perturbations in temperature and vor- 
ticity to break the axisymmetry of the initial state. Each 
simulation was integrated for at least 40,000 Earth days. 
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In the statistically steady states, the global-mean outgo- 
ing thermal radiative flux is within 1% of the sum of the 
global-mean solar radiative flux and the imposed intrin- 
sic heat flux. The vertically integrated Rayleigh drag on 
the zonal flow approximately balances the vertically in- 
tegrated total (mean plus eddy) angular momentum flux 
convergence. The circulation statistics shown are com- 
puted from flow fields sampled 4 times daily and av- 
eraged over 1500 days. They were first computed on 
the GCM's sigma surfaces, with the appropriate sur- 
face pressure-weighting of the averages (Schneider and 
Walker 2006), and then interpolated to pressure surfaces 
for display purposes. 
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